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HIGHER  HARMONIC  GENERATION  IN  THE 
INDUCED  RESONANCE  ELECTRON  CYCLOTRON  MASER 


I .  INTRODUCTION 

Generation  of  intense  radiation  in  the  microvave  regime  through 

electron  cyclotron  interaction  vas  proposed  independently  by  a  number  of 
1-4 

researchers  in  the  late  1950s.  Electrons  gyrating  in  resonance  with  the 

radiation  field  experience  a  bunching  in  the  relative  wave  particle  phase 

through  the  dependence  of  the  cyclotron  frequency  on  the  relativistic  mass. 

High  amplification  of  the  radiation  field,  known  as  masing  action,  results 

for  radiation  frequencies  slightly  above  the  electron  cyclotron  frequency. 

5- 18 

Electron  cyclotron  masers  ,  also  called  gyrotrons,  have  demonstrated 
efficient  high  power  generation  capability  at  the  centimeter  wavelengths. 

Electron  cyclotron  instabilities  also  occur  in  ionospheric  and  astrophysic 

,  19,20 

plasmas 

A  variety  of  potential  applications,  such  as  advanced  accelerators, 

heating  of  fusion  plasmas,  short  wavelength  radar  and  spectroscopy,  call 

for  generation  of  intense  radiation  at  even  shorter  wavelengths  in  the 

millimeter  and  submillimeter  range.  In  a  closed  resonator,  the  shortest 

wavelength  for  single  mode  operation  is  tied  to  the  transverse  dimension  of 

the  cavity.  Operation  at  radiation  wavelengths  shorter  than  the  transverse 

21 

dimensions  will  result  in  a  multimode  excitation  due  to  the  small 


frequency  separation  among  cavity  eigenmodes.  This  limitation  in  the 

22  23 

wavelength  is  considerably  relaxed  in  the  quasi-optical  maser  ’ 
operating  in  an  open  resonator  that  offers  much  improved  frequency 


separation. 

Considerable  attention  has  been  given  lately  to  the  operation  at 

Doppler  upshifted  f requencies"^-'^  resulting  from  a  finite  wave  number 

in  the  direction  of  the  electron  beam  propagation.  The  operation  frequency 

w,  defined  by  the  resonance  condition  oo-kz-Q  =0,  is  given  bv 
J  z  c 
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•  =  2C(1  '  nz02)  ’ 

2  -1/2 

where  y  is  the  relativistic  factor  y  =  (1-0  )  ,  0  =  v/c,  =  kzc/w  is 

the  parallel  index  of  refraction  and  S  =  3  /y  is  the  relativistic 

co 

cyclotron  frequency  with  2q  =  eBQ/mc.  For  nz  =  0^  =  1  the  frequency  is 
2 

boosted  to  2yz  times  the  electron  cyclotron  frequency  with 
2  -1/2 

Yz  =  (l-@z  )  •  Note  that  both  parallel  and  perpendicular  kinetic  energy 

of  the  electrons  feed  the  instability  in  case  of  a  tilted  resonator.  So 

far  plane  wave  configurations  in  simple  geometry,  also  referred  to  as  the 

24  24-27 

cyclotron  autoresonance  maser  (CARM),  have  been  analyzed  in 

conjunction  with  Doppler  upshifting  of  the  radiation  frequency. 

28-31 

The  induced  resonance  electron  cyclotron  (IREC)  maser  ,  shown  in 
Fig.  1,  operates  at  Doppler  upshifted  frequencies,  utilizing  at  the  same 
time  the  advantages  of  the  open  resonators.  Each  resonator  forms  an  angle 
a  with  the  direction  of  the  electron  beam  along  the  external  magnetic 
field.  The  index  of  refraction  nz  =  cosa  is  adjustable  by  varying  the  angle 


between  the  resonators  and  can  be  chosen  to  minimize  the  effects  of  finite 

beam  thermal  spreads.  For  operation  at  the  optimum  refraction  index  the 

efficiency  is  relatively  insensitive  to  the  beam  energy  spread.  The 

sensitivity  to  the  beam  pitch  angle  spread  can  also  be  minimized.  The 

interaction  length  inside  the  resonator  that  maximizes  efficiency  is  of  the 

14 

order  of  half  a  bounce  distance  for  the  trapped  particles. 

As  the  available  magnetic  field  limits  the  maximum  operation  frequency 

31  32 

for  given  y,  operation  at  higher  harmonics  *  becomes  very  attractive. 

The  magnetic  field  required  to  produce  radiation  at  a  given  frequency  is 


reduced  to  a  fraction  1/(N  +  1)  for  operation  at  the  Nth  harmonic. 

33 , 34 

Operation  at  higher  harmonics  has  been  analyzed  for  the  conventional 
35  36 

and  quasioptical  ’  gyrotron,  showing  that  considerable  efficiency  can  be 


achieved  at  the  first  few  harmonics.  Previous  IREC  studies  have  considered 
operation  at  the  fundamental  frequency  and  small  Larmor  radius,  relevant  to 
the  case  of  small  resonator  angle  and  small  pitch  angle  0  =  v1/vz-  In  the 
present  work  we  analyze  operation  at  any  given  harmonic  N  and  arbitrary 
resonator  angle  including  finite  Larmor  radius  effects.  The  effects  of  the 
Gaussian  radiation  profile  are  retained  as  well.  A  set  of  slow  time  scale 
equations  of  motion  is  derived  by  averaging  over  the  cyclotron  period  time 
scale.  The  small  signal  efficiency  is  determined  analytically.  Nonlinear 
efficiency  is  determined  by  numerical  integration.  It  is  found  that  the 
efficiency  for  the  first  few  harmonics  is  comparable  to  that  for  the 
fundamental  in  the  same  parameter  regime.  This  is  feasible  because 
saturation  occurs  at  larger  radiation  amplitude  with  higher  harmonics.  The 
effects  of  energy,  pitch  angle  and  guiding  center  spread  are  also  studied. 
An  optimum  resonator  angle  a  exists  for  a  given  set  of  parameters 
minimizing  the  effects  of  finite  beam  thermal  spreads.  Efficiency 

enhancement  can  be  achieved  by  properly  tapering  the  external  magnetic 

29 

field  ,  inducing  an  extended  wave  particle  resonance. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Section  II, 
we  describe  the  field  in  the  resonator  and  we  obtain  the  fast-time  averaged 
equation  of  motion.  In  Section  III,  the  small  signal  efficiency  and  the 
start-up  current  required  to  trigger  the  oscillations  in  the  resonator  are 
calculated.  In  Section  IV,  we  discuss  briefly  the  effects  caused  by  the 
finite  thermal  spreads  in  the  electron  energy  and  pitch  angle.  In  Section 
V,  we  integrate  numerically  the  equations  of  motion  using  velocity 
distributions  with  finite  energy  and  pitch  angle  spread  as  veil  as  guiding 
center  distribution  in  the  transverse  direction.  The  nonlinear  efficiency 
is  computed  for  the  first  few  harmonics. 


II.  FIELD  MODELING  AND  PARTICLE  DYNAMICS 
The  configuration  for  the  induced  resonance  electron  cyclotron  (IREC) 
maser  is  shown  schematically  in  Fig.  1.  The  interaction  cavity  is  formed 
by  the  two  quasi-optical  resonators  intersecting  at  an  angle  2a  where  a  is 
the  angle  relative  to  the  external  magnetic  field  Bq  along  the  z-axis.  The 
electron  beam  also  propagates  along  z.  The  total  vector  field  is  the 
superposition  of  the  two  resonant  fields 


A(x',y',z';t  )  +  A(x' ,y' ,z' ; t)  ,  (1) 
ot  -ot 

where  A(x' ,y' , z' ; t )  are  eigenmodes  of  the  Fabry-Perot  type  resonator.  Here 
we  consider  the  lowest  order  Gaussian  TEMqq  modes  linearly  polarized  along 
the  y-axis 


2  .2 


A( 


zoo  .  <v  z  V^htl  •  u<z'>) 


exp  (i(kz'-<ot)]  +  exp  [-i(kz'  +  ut')]>  +  c.c., 


(2) 


where 


1  1 
q(z' )  =  R(z' ) 


itw^(z' ) 


R(z') 


.2,1/2 


w(z')  -  vo[l  ♦  (f^-)  ]  ,  u(z')  =  tan^g-), 


w  is  the  beam  waist  at  the  center  of  the  resonator,  A  is  the  wavelength 
o 

2 

and  Zq  =  mWA  t^ie  Rayleigh  length. 

The  coordinates  (x',y',z')  have  the  z'-axis  aligned  with  each 

£01 

resonator  and  are  related  to  (x,y,z)  by 


'x'j 

'coso  ±sina^ 

.z'  ; 

,+  sina  cosa/ 

Zj 

y'  =  y  - 


4 


The  radiation  field  is  assumed  at  temporarily  steady  state.  The 
Gaussian  width  wq  for  the  radiation  envelope  is  much  larger  than  the 
radiation  wavelength  A  and  the  beam  spot  size  b.  The  Rayleigh  length  Zq  is 
typically  much  longer  than  the  interaction  length  L  -  vQ/sina  in  the 
z-direction.  In  the  vicinity  of  the  beam  x  ~  y  -  b,  and  within  the 

o 

interaction  regime  |z|  -  L,  we  have  b/L  -  z/Zq  ~  e  «  1,  b/ZQ  -  t. 

2 

Combining  Eqs.  (1)  -  (3),  and  dropping  terms  of  order  e  ,  the  forward 
propagating  component  of  the  total  field  near  the  interaction  area  is 
expressed  as 

2  .  2 
z  sin  a 

2 

At(x,y,z;t)  =  eyA0e  w0  cos  k^x  C0S  (kzz-ttt),  (4) 

where  k^  =  (w/c)  sin  a,  k  =  (ai/c)  cos  at. 

In  short,  Gaussian  effects  in  the  transverse  direction  of  order 
2  2  2 

b/L  ~  z  have  been  omitted  while  the  electrons  experience  a  Gaussian 
envelope  of  effective  width  L  =  wo/sina  in  the  z-direction.  Only  the 
forward  propagating  wave  phase  is  considered  in  Eq.  (4)  since  the 
synchronous  interaction  of  an  electron  with  the  backward  component  occurs 
at  down-shifted  cyclotron  frequency,  of  small  practical  interest. 

We  use  the  guiding  center  description  for  the  particle  orbits 


5 


time  scale,  the  slow  time  scale  nonlinear  relativistic  equations  of  motion 
are  cast  in  the  form 

dx  ^ 

di*  =  ‘  2^  i~  a  JN(klp)  sin  *N  cos  gN’  (6a) 

o  z 


dZs_ 

dz 

Ckl  U1 

—  Y  ~  a 

0  z 

J^,(klp)  cos  v|»N  cos  gN, 

(6b) 

du^ 

dz 

.  ic  ) 

lcuz  zj 

a  J^,(kiP)  sin  *N  sin  gN> 

(6c) 

du 
_ z 

dz 

U1 

kz  r  a  JN(klp)  sin  *N  sin  *N- 
z 

(6d) 

dz 

vto  N2 

-k*' h 

fyw  .  1  JN(klp) 

l^r-kzja  k,p  cos^smgN 

Z  J. 

kl 

-  N  —  a  JN(klp)  cos  %  sin  gN- 
z 

(6e) 

In  Eqs.  (6a)  to  (6e)  the  prime  (')  signifies  the  Bessel  function 

derivative  in  respect  to  the  argument,  u  is  the  normalized  momentum 

2  2  1/2 

u  ■  p/mc  *  yv/c,  y  is  the  relativistic  factor  y  =»  (1  +  u^  *  u^)  , 

^  =  k^z  -  wt  +  NC  +  Nn/2  is  the  relative  phase  between  the  field  and  the 

particle,  gN  =  k^  xg  -  Nn/2  carries  the  dependence  on  the  guiding  center 

2  2 

position,  a(z)  =  aQ  exp(-z  /L  )  is  the  normalized  radiation  amplitude  with 
2 

aQ  =  eAQ/mc  and  the  detuning  in  frequency  Aw  is  given  by 

Aw  =  w(l-n  3  )  -  NS2  /y.  (7) 

zz '  o 

It  is  the  dependence  of  the  detuning  Aw  on  the  particle  energy  through  the 
relativistic  correction  y  that  causes  the  phase  bunching  and  the  radiation 


6 


amplification  in  the  linear  regime.  The  evolution  of  y  is  found  combining 
Eqs.  (6c)  to  (6e), 


I 


ir 

dz 


u. 

c  r  3  JN(klp)  Sin  *N  Sin  gN- 
z 


(8) 


In  performing  the  fast  time  averaging  to  obtain  Eqs.  (6)-(8)  we  have 

assumed  that  the  particles  always  remain  close  to  resonance  with  a  single 

harmonic  N,  i.e.,  <o(l-nz0z)  -  N2Q/y  =  0.  Therefore  the  change  in  energy  Ay 

and  consequently  the  radiation  amplitude  aQ  cannot  exceed  a  certain  limit. 

When  aQ  is  very  large,  the  particles  may  also  experience  resonant  effects 

from  nearby  harmonics  N  ±  1.  The  validity  conditions  for  a  single  resonant 

harmonic  are  satisfied  in  the  parameter  regime  under  consideration. 

The  nonlinear  system  of  differential  equations  ( 6 )  —  < 8 )  cannot  be 

solved  analytically  in  terms  of  elementary  functions,  except  in  some 
37  38 

special  cases  ’  .  We  resort  to  numerical  integration  in  order  to  examine 


the  nonlinear  behavior  while  the  small  signal  analysis  is  done  by 
perturbation  theory. 


III.  SMALL  SIGNAL  EFFICIENCY 


One  of  the  issues  concerning  intense  microwave  generation  is  the 
intrinsic  efficiency  h  of  the  interaction,  defined  as 


h 


<Yf-Yo> 
<r0-1  > 


<Yo-1>"1  JdPJ/(P0)  J  dgQ  J  d*Q  |  dz  g, 


(9) 


where  <  >  signifies  the  average  over  the  initial  distribution  in  phase 

space,  and  y  is  a  function  of  the  initial  conditions  r(z»  p .  .  p  ,  ♦  ,  g  ) 

with  y  and  gQ  »  g^(-®).  We  compute  the  small  signal  efficiency 

in  order  to  determine  the  beam  current  required  to  overcome  losses  and 

start  the  cavity  oscillations.  After  obtaining  the  linearized  solutions  of 

Eqs.  (6)  to  (8),  we  substitute  them  into  the  integrant  in  the  right-hand 

side  of  Eq.  (9).  The  evaluation  of  the  final  result  is  considerably 

simplified  by  taking  the  phase  space  average  over  before  the  spatial 
22  36 

integration  ’  over  z.  Expanding  the  products  of  trigonometric  terms 
inside  the  integrant  into  sums,  averaging  over  <pQ  and  extending  the  limits 
of  the  z-integration  to  ±«,  we  obtain  the  linear  efficiency  in  terms  of  the 
initial  conditions 


2  2 

a  £  .  .2  /  i  j  too"  \ 

n  *  4  r0770-i~  (jn<VJ  exp  (-  I  *  IT  J  <sin  *0> 


•  f  N‘J  (s  ) 

n  0  -1  1  +  — 777 

I  z  zo  )[  soJ^(so) 


Vn(so^ 


+  JN<so>  ) 


+  90 
0  zo 


NW  1 

eonz  +  J^(so)  nl 


(  n  n  n  <COS  g  >  SJ.Js)) 

-2«  2,,  2v  °o  o  N  o' 


.  .  2  .  J'(s  ) 

<sin  g  >  No 
so 


.Aw  N  ~  n  ,Aw  .2) 

f— )  -  n  0  02£2(— )  1,  (10) 

l  W  J  Z  ZO  0  {  w  J  I 


where  sq  *  k^,  gQ  =  k±xgQ  -  Nn/2,  0iQ  =  v^/c,  0ZQ  =  v^/c,  9o  =  v^/v^ 
and  <f>  =  (l/2n)  f2rt  d(k^x^)/  is  the  average  over  the  initial  guiding 
center  position. 


Ve  have  chosen  to  express  q  in  terms  of  the  parameters  Awq/w  and 
£  =  wx,  where  T  =  yQL/cu^Q  is  the  transit  time  through  the  interaction 
regime.  The  argument  £Awq/w  in  the  exponential  is  equal  to  Awqt,  the 
linear  advance  in  the  relative  phase  A^  over  the  interaction  regime. 

The  sign  is  determined  by  the  angular  bracket  on  the  right-hand  side  of 
Eq.  (10).  Treating  the  bracket  as  a  quadratic  form  in  Awq/w  and  keeping 
the  lowest  order  contribution  in  k^p,  we  find  that  the  rejime  for  positive 
efficiency  is  given  approximately  by 


(N+l )  fl 

\ 

'i  2  1 

L  -  n  0  -  9  n  0 

z  zoj  o  z^zo. 

Aw 

,  0  ^ 

2  2 
eio(1  -  V 

] 

t. 

2  2 

L  -  n^  0,  r 
z  lo  0 

-  ^  ^ 

U) 

n  0  92 

z  zo  0 

The  upper  limit  in  Aw  /w  is  due  to  the  finite  n  and  results  from  the 

o  z 

2 

negative  contribution  of  the  quadratic  term  (Aw/w)  that  overtakes  the 
positive  contribution  of  the  linear  terms  for  small  angle, 


si 


n2<x 


n  8  02 

z  zo  o 


For  typical  operation  parameters  we  have  t  >>  1  and  Awq/w  <<  1.  In 

order  to  estimate  the  maximum  efficiency  within  the  positive  regime,  we 

parametrize  Eq.  (10)  as  a  function  of  C  =  £Awq/w  and  look  for  the  zeros  of 

the  cubic  equation  resulting  from  dh/dC  =  0.  If  the  angle  a  is  not  too 
2  2 

small,  sin  a  >>  l/£0,  where  £,  >>  1,  then  the  maximum  occurs  at  C  =  1, 
yielding 


.2  P3 


max 


it 

8 


kv,')' 


lo 


(y!) 


y  sin 
o 


exp 


(-  \), 


(12) 


2  2  2  2 

where  £  =£(l-n)=(wy  w/cu  )  is  independent  of  a.  The  small 

o  z  oo  zo  K 

signal  efficiency  increases  with  decreasing  angle  a  (increasing  index  of 
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refraction  nz)  provided  that  sina  satisfies  the  inequality  above  Eq.  (12). 
When  sina  is  too  small,  an  exact  solution  of  the  cubic  equation  for  £  is 
required  and  Eq.  (12)  is  invalid. 

We  now  calculate  the  start-up  beam  current  using  the  small  signal 
efficiency.  Amplification  of  the  electromagnetic  radiation  is  possible  if 


s  dE 

nPb  >  dt  =  Q  E’ 


(13) 


2  2  2  2  2  b  2 

where  E  »  (l/4)w^  L^a^fw  /c  )  (me  /|e|  ;  is  the  total  energy  stored  inside 
both  resonators,  is  the  resonator  length,  dE/dt  is  the  combined 
refraction,  diffraction  and  transmission  losses  and  Q  is  the  quality  factor 
for  the  cavity.  Combining  Eqs.  (12)  and  (13)  we  obtain 


Y_ ( Y  -1 )  sin  a  exp  (-  \)  25  L_ 

j  y  y  2  20  o  o _  _  Z  me  R 

s  b  a  2 

0lo 


11  [  Pq  )  1 


I  i2  v  Q' 

le  I  o 


(H) 


For  small  Larmor  radius  k,  p  «  1  the  start-up  current  I  increases  very 
quickly  with  the  harmonic  N,  I  «  22N[ <N-1) ! 1 2/ ( k± pq ) ~ 1 ^ .  It  is 
therefore  desirable  to  operate  at  k^pQ  >  1  in  order  to  have  good  coupling 
to  the  resonator  modes  and  low  start-up  currents.  In  this  case  the  start¬ 
up  current  scales  roughly  as  the  inverse  maximum  of  the  Bessel  function 
1/3 

derivative,  I  *  N  ,  increasing  mildly  with  the  harmonic  N.  We  can 
achieve  harmonic  selection  by  choosing  k^pQ  near  a  maximum  of  JfJ'(k^po)  for 
the  intended  Nth  harmonic. 

Expression  (14)  for  the  start-up  current  was  derived  using  coherent 
resonator  modes.  These  modes  have  evolved  from  an  initial  noise  background 
of  spontaneous  cyclotron  radiation.  During  the  spontaneous  emission  stage 
preceding  coherency  most  of  the  emitted  radiation  is  contained  within  a 
cone  of  angle  1/y  around  the  direction  of  the  particle  velocity  v.  Since  v 
makes  an  angle  9  =  tanfv^/v^)  with  the  magnetic  field  the  condition 
|0  -  a|  <  1/y  must  be  met  to  avoid  excessive  losses  during  the  start-up 


IV.  THERMAL  EFFECTS 


One  of  the  important  features  associated  with  the  IREC  maser  is 
choosing  the  index  of  refraction  appropriately  to  minimize  the  effects  of 
the  electron  beam  thermal  spreads.  Spreads  in  the  initial  electron 
momentum  will  cause  a  spread  in  the  detuning  parameters  dwo  among  different 
particles.  This  in  turn  will  cause  an  accelerated  mixing  in  phase  space 
opposing  the  nonlinear  phase  bunching  and  reducing  in  efficiency.  According 
to  Eq.  (7)  a  standard  deviation 


.2  „l/2 


«*•><,  ■  [(§r)  i6l]  -  [<-v  nVo0*o>X*<[,v,'V25'jJ] 


1/2 


(15) 


in  the  initial  detuning  results  from  a  beam  distribution  with  velocity 
deviations  S3Z  and  53^.  The  choice  of  resonator  angle 


N2  y  0 
O  O  20 

n  =»  cosa  =  - 

z  co 


(16) 


causes  the  minimum  initial  spread  in  for  any  beam  thermal  spreads.  The 
minimum  spread  from  Eqs.  (15)  and  (16)  can  also  be  expressed  in  terms  of 
the  pitch  angle  spread  5 0q  and  the  energy  spread  5yQ.  Then,  the 
requirement  for  small  phase  mixing  among  various  particles  over  the 
interaction  length  L,  namely  S(Aw)oL/c3z  <<  rc,  suggests 


50 
_ o 

9 
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where  Nc  =  QqL/ (2tiyoc320)  is  the  approximate  number  of  cyclotron  gyrations 
within  the  interaction  length.  The  beam  thermal  spread  requirements  become 
more  stringent  with  increasing  harmonic  N. 

In  the  nonlinear  operation  regime  the  electrons  get  trapped  in  the 
wave  potential^’  and  execute  synchrotron  oscillations  in  phase  space  in 
a  similar  manner  as  in  the  conventional  gyrotrons.  The  phase  mixing  among 


various  particles  is  now  determined  by  the  dependence  of  the  trapped 
particle  synchrotron  period  on  the  various  parameters.  The  efficiency 
deterioration  involves  more  factors  than  the  initial  spread  in  the  detuning 
which  is  the  dominant  source  of  phase  mixing  only  in  the  small  signal 
regime  aQ  <<  1.  Analytic  predictions  similar  to  Eq.  (16)  are  hard  to  make 
in  the  nonlinear  case.  Our  numerical  results  show  that  when  the  index  of 
refraction  is  optimized  according  to  (16)  the  nonlinear  efficiency  is 
practically  insensitive  to  the  energy  spread  5y  /y  • 


V.  NUMERICAL  RESULTS 


In  this  section,  we  investigate  various  aspects  of  the  nonlinear 

performance  by  numerically  integrating  Eqs.  (6a)  to  (6e).  We  consider  an 

electron  beam  of  r  =  2.5  (~  0.75  MeV)  with  0.  =  1/y  and  0  =  0.825  in 

o'  io  o  20 

a  guide  magnetic  field  of  strength  Bo  =  40  kG.  The  appropriate  index  of 

refraction  to  minimize  the  effect  of  energy  spread  is,  according  to  Eq. 

(16),  n  »  0.982  which  corresponds  to  an  angle  a  =  11°.  The  frequency  is 

upshifted  by  a  factor  N/(l  -  n02Q)  =  5.26  N  times  the  relativistic 

cyclotron  frequency,  and  corresponds  to  a  wavelength  of  0.41  mm  for  the 

third  harmonic  N  =  3  and  0.31  mm  for  the  fourth  harmonic  N  =  4.  The 

radiation  spot  size  wq  is  0.50  cm  and  the  Rayleigh  lengths  are  14.3  cm  and 

25.8  cm  for  N  »  3  and  N  »  4  respectively.  We  consider  a  uniform  guiding 

center  distribution  in  the  interval  0  <  k.x  <  2it. 

1  g 

Curves  of  efficiency  n  versus  aQ  for  various  values  of  the  detuning 

AdWw  are  plotted  for  N  =  3  and  N  =  4  in  Figs.  2  and  3  respectively.  These 

results  correspond  to  a  cold  beam  without  thermal  spreads.  We  find  the 

efficiency  for  the  first  few  harmonic  comparable  to  the  efficiency  for  the 
29 

fundamental  under  the  same  operation  parameters.  Efficiency  saturation 
occurs  for  larger  amplitude  aQ  compared  to  the  operation  at  the 
fundamental.  Figure  4  shows  the  effects  of  finite  beam  quality  on 
efficiency  when  the  electron  beam  has  either  a  spread  in  the  pitch  angle  or 
a  spread  in  energy.  We  plot  the  ratio  of  the  thermalized  efficiency  h  over 
the  cold  beam  efficiency  h  for  the  third  harmonic  N  =  3  at  fixed  amplitude 
aQ  =  0.20.  Curve  (a)  for  zero  energy  spread,  Ayo/yq  =  0,  shows  that  the 
half  width  in  the  pitch  angle  spread  that  reduces  efficiency  by  50%,  is 
equal  to  A9q/9o  =  2%.  Curve  (b)  for  zero  pitch  angle  spread,  A9o/9o  =  0, 
shows  that  the  half  width  in  energy  spread  is  A yq/yo  =  13%.  Efficiency 
tends  to  be  more  sensitive  on  the  spread  in  the  pitch  angle  than  the  spread 


in  energies;  therefore,  we  may  simulate  thermal  effects  by  including  only 
pitch  angle  spreads,  cutting  down  on  computing  time. 

Given  that  the  large  signal  efficiency  depends  on  few  parameters, 
predictions  about  optimum  operation  at  maximum  efficiency  are  hard  to  make. 
One  anticipates  maximum  efficiency  when  the  transit  time  through  the 
interaction  regime  is  about  equal  to  half  the  synchrotron  period  for  a 
trapped  particle.  An  optimum  interaction  length  in  the  z  direction 
Lz  =  2L  =  2wQ/sina  corresponds  to  a  given  synchrotron  period  o^,  which,  in 
turn,  depends  on  the  five  parameters  aQ,  rQ,  9Q,  and  cosa.  This  is 
illustrated  in  Fig.  5,  which  shows  efficiency  as  a  function  of  the  traveled 
distance  z  for  three  different  Gaussian  profiles  corresponding  to  different 
radiation  spot  sizes  wq,  keeping  the  other  parameters  fixed.  In  curve  (a) 
the  interaction  length  is  less  than  half  the  bounce  distance  =  nc0  /w. 

and  the  electrons  exit  the  resonator  before  reaching  the  point  of  lowest 

energy  in  their  trajectories.  In  curve  (b)  we  have  a  good  matching  of  Lz 
with  achieving  the  highest  efficiency.  In  curve  (c)  Lz  is  larger  than 
and  the  electrons  overshoot  the  point  of  minimum  energy,  gaining  energy 
back  from  the  wave  and  reducing  efficiency. 

From  the  practical  point  of  view,  one  would  like  to  optimize  the 

design  parameters  of  the  resonator  a  and  L  for  a  given  beam  energy  yq  and 

2 

pitch  angle  9q  under  the  maximum  energy  load  aQ  sustained  by  the  cavity. 

We  already  picked  the  operation  angle  cosa  so  as  to  minimize  the  effects 

the  beam  energy  spread.  In  Fig.  6  we  show  the  efficiency  as  a  function  of 

the  interaction  length  L  by  varying  the  spot  size  w  and  keeping  all  other 

parameters  fixed.  The  upper  curve  shows  the  nonlinear  efficiency  for  a 

monoenergetic  electron  beam  of  infinitesimal  spot  size  b  <<  w  ,  and  a 

uniform  spread  in  the  initial  phase  0  <  4^  <  2 n.  A  uniform  guiding  center 

distribution  0  <  k,x  <  2 ri  is  included  in  the  second  curve.  The  resulting 
-  I  g  -  s 
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efficiency  reduction  is  no  more  than  30 7.  indicating  that  some  bunching  in 
the  guiding  center  position  takes  place  as  well.  The  addition  of  a  27. 
energy  spread  &ro/Y0  with  zero  pitch  angle  spread  does  not  reduce 
efficiency  considerably  in  the  fourth  curve.  A  27.  spread  in  the  pitch 
angle  50q/9o  with  zero  energy  spread  has  a  more  pronounced  effect  on 
efficiency  shown  in  the  lowest  curve  (d).  The  overall  picture  shows  that, 
for  the  parameters  chosen,  efficiency  has  a  weak  dependence  on  the 
interaction  length  L  falling  off  slowly  after  an  optimum  length  of 
L  -  4  cm. 
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Figure  1.  Schematic  illustration  of  the  Induced  Resonance  Electron 


Cyclotron  Maser. 


Plots  of  efficiency  h  versus  the  normalized  radiation 
amplitude  aQ  at  the  third  harmonic  N  =  3  and  for  detuning 
parameters  AaWw  equal  to  (a).  0.025  (b).  0.050  (c).  0.075  and 
(d).  0.100  respectively.  A  cold  beam  of  uniform  guiding 
center  spread  is  considered.  The  simulation  parameters  are 


Figure  3.  Plots  of  efficiency  h  versus  the  normalized  radiation 
amplitude  aQ  at  the  fourth  harmonic  N  =  4  for  detuning 
parameter  u  equal  to  (a).  0.025  (b)  0.037  and  (c)  0.050 

respectively.  The  other  parameters  are  the  same  as  in  Fig.  2 


gure  4. 
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Dependence  of  the  efficiency  on  beam  thermal  spreads.  Shown 

is  the  ratio  of  thermalized  to  cold  beam  efficiency  n/hQ  as  a 

function  of  (a)  pitch  angle  spread  S9q/9o  with  =  0  and  (b) 

energy  spread  5yq/yo  yith  S9q  =  0.  The  parameters  are  the 

same  as  in  Fig.  2  with  Aw  /w  =  0.075. 
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Figure  5.  Plots  of  efficiency  n  versus  travelled  distance  z  inside  the 

resonator  for  various  interaction  lengths  corresponding  to 

different  radiation  spot  sizes  v  .  The  center  of  the 

o 

resonator  is  at  z  =  0.  is  equal  to  (a)  2.64  cm  (b)  4.24  c: 

and  (c)  7.41  cm.  The  parameters  are  the  same  as  in  Fig.  2 

vith  A u  /<*>  =  0.075. 
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